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Abstract 

This paper is a review of previous works on the stochastic background 
of gravitational waves (SBGWs) which has been discussed in various peer- 
reviewed journals and international conferences. The SBGWs is analyzed 
with the aid of the Wilkinson Microwave Anisotropy Probe (WMAP) 
data. We emphasize that, in general, in previous works in the literature 
about the SBGWs, old Cosmic Background Explorer (COBE) data were 
used. After this, we want to face the problem of how the SBGWs and 
f(R) gravity (where f(R) is a function of the Ricci scalar R) can be related, 
showing, vice versa, that a revealed SBGWs could be a powerful probe 
for a given theory of gravity. In this way, it will also be shown that the 
conform treatment of SBGWs can be used to parametrize in a natural 
way f(R) theories. Some interesting examples which have been recently 
discussed in the literature will be also analysed. 

The presence and the potential detection of the SBGWs is quite impor- 
tant in the framework of the debate on high-frequency gravitational waves 
(HFGWs) too. Recently, the importance of HFGWs has been emphasized 
in some papers in the literature. 

PACS numbers: 04.80.Nn, 04.30.Nk, 04.50.+h 

Keywords: gravitational waves; extended theories of gravity; stochastic back- 
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1 Introduction 



The accelerated expansion of the Universe, which is today observed, shows that 
cosmological dynamic is dominated by the so called Dark Energy which gives 
a large negative pressure. This is the standard picture, in which such new 
ingredient is considered as a source of the right hand side of the well known 
Friedman-Robertson- Walker field equations [T] . It should be some form of non- 
clustered and non-zero vacuum energy which, together with the clustered Dark 
Matter, drives the global dynamics. This is the so called "concordance model" 
(Lambda Cold Dark Matter, ACDM) which gives, in agreement with the Cos- 
mic Microwave Background Radiation (CMBR), Large Structure Scale (LSS) 
and Supernovae la (SNela) data, a good framework of the today observed Uni- 
verse, but presents several shortcomings as the well known "coincidence" and 
"cosmological constant" problems p] . 

An alternative approach is changing the left hand side of the field equations, 
seeing if observed cosmic dynamics can be achieved extending general relativity 
[HOHHS]. In this different context, it is not required to find out candidates for 
Dark Energy and Dark Matter, that, till now, have not been found, but only 
the "observed" ingredients, which are curvature and baryon matter, have to be 
taken into account. Considering this point of view, one can think that gravity 
is not scale- invariant and a room for alternative theories is present [3 [8] . 

In this picture, even the sensitive detectors for gravitational waves, like bars 
and interferometers (i.e. those which are currently in operation and the ones 
which are in a phase of planning and proposal stages) [HI El [HI [HI 031 HH 
1 151 116j . could, in principle, be important to confirm or to rule out the physical 
consistency of general relativity or of any other theory of gravitation. This 
is because, in the context of Extended Theories of Gravity, some differences 
between General Relativity and the others theories can be pointed out starting 
from the linearized theory of gravity [SJ [T71 [T51 [TO1 [52] ■ 

In principle, the most popular Dark Energy and Dark Matter models can 
be achieved considering f(R) theories of gravity [2][3j|4], where R is the Ricci 
curvature scalar. In the tapestry of Extended Theories of gravity, it is quite 
important to understand why we are going to select f(R) theories of gravity 
among several existing modified theories of gravity |53j . The motivation on 
such a choice have been recently well emphasized in the interesting analysis 
in [54]. Thus, following the review [54], we recall that some f(R) theories of 
gravity are not excluded by requirements of Cosmology and Solar System tests, 
i.e are viable. On the same time, they have quite rich cosmological structure: 
they may naturally describe the effective (cosmological constant, quintessence 
or phantom) late-time era with a possible transition from deceleration to ac- 
celeration thanks to gravitational terms which increase with scalar curvature 
decrease. In other words, f(R) theories are well tuned with an accelerating 
standard FRW universe. These theories may be also rewritten as scalar-tensor 
gravity (Einstein frame) with the aid of conform transformations. Starting by 
these considerations, in [M] two specific viable models have been considered: the 
one with 1/R and R m terms, and another with lni? terms. It has been shown 



2 



that both of the models may lead to the (cosmological constant or quintessence) 
acceleration of the universe as well as an early time inflation. Moreover, the first 
model seems to pass the Solar System tests, i.e. it has the acceptable newto- 
nian limit, no instabilities and no Brans-Dickc problem (decoupling of scalar) 
in scalar-tensor version [S3]. These versions of f(R) theories may also describe 
cosmic speed-up at late universe. 

f{R) theories can be taken into account also in the analysis of the SBGWs 
which, together with the CMBR, would carry, if detected, a huge amount of 
information on the early stages of the Universe evolution [201 [HI E2])- A key 
role for the production and the detection of this graviton background is played 
by the adopted theory of gravity [3TJ [S3] in this case too. 

In the second section of this review, the SBGWs is analysed with the aid of 
the WMAP data [20[ [22J [23] . We emphasize that, in general, in previous works 
in the literature about the SBGWs, old COBE data were used (see [24 l l25 l l26 l 
EH EH [29] for example). 

In the third section we want to face the problem of how the SBGWs and 
f{R) gravity can be related, showing, vice versa, that a revealed SBGWs could 
be a powerful probe for a given theory of gravity. In this way, it will also be 
shown that the conform treatment of SBGWs can be used to parametrize in 
a natural way f(R) theories. In this framework, we will discuss like particular 
case [S3] some interesting viable f(R) theories that have been analysed in refs. 
[SJ (SSJ [SB] . In this way, a discussion which started in [S7] will be integrated. 

The presence and the potential detection of the SBGWs is quite important 
in the framework of the debate on high-frequency gravitational waves (HFGWs) 
too. Recently, the importance of HFGWs has been emphasized in some papers 
in the literature gEl SSI EOj [ST] . 

We emphasize that this is a review of previous works on the SBGWs which 
have been discussed in various peer-reviewed journals and international confer- 
ences. 

2 Constrains on the stochastic background of 
gravitational waves from the WMAP data 

From our analysis, it will result that the WMAP bounds on the energy spectrum 
and on the characteristic amplitude of the SBGWs are greater than the COBE 
ones, but they are also far below frequencies of the earth-based antennas band. 
At the end of this section a lower bound for the integration time of a potential 
detection with advanced LIGO is released and compared with the previous one 
arising from the old COBE data. Even if the new lower bound is minor than 
the previous one, it results very long, thus for a possible detection we hope in 
the Laser Interferometric Space Antenna (LISA) and in a further growth in the 
sensitivity of advanced projects. 

The strongest constraint on the spectrum of the relic SBGWs in the fre- 
quency range of ground based antennas like bars and interferometers, which is 
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the range 10Hz < f < 10 4 Hz, comes from the high isotropy observed in the 
CMBR. 

The fluctuation AT of the temperature of CMBR from its mean value To = 
2.728 K varies from point to point in the sky [501 HH E3] , and, since the sky can 
be considered the surface of a sphere, the fitting of AT is performed in terms 
of a Laplace series of spherical harmonics 

AT 00 1 

— (f)) = E 0«milm(A) (1) 

I— 1 m— — l 

and the fluctuations are assumed statistically independent (< ai m >= 0, 
< ai m a* Vm , >= Ci5u>5 mm >). f) denotes a point on the 2-sphere in Eq. Q}, while 
the ai m are the multipole moments of the CMBR. For details about the defini- 
tion of statistically independent fluctuations in the context of the temperature 
fluctuations of CMBR see [H US]. 

The WMAP data [HI [53] permit a more precise determination of the rms 
quadrupole moment of the fluctuations than the COBE data 
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Q rms = T(J2 = 8 ± ^K, (2) 

m=-2 

while in the COBE data we had [M [M [271 



Q rms = 1 i..r;;n//A. (3) 

A connection between the fluctuation of the temperature of the CMBR and the 
SBGWs derives from the Sachs- Wolfe effect [101 ISO]- Sachs and Wolfe showed 
that variations in the density of cosmological fluid and GWs perturbations re- 
sult in the fluctuation of the temperature of the CMBR, even if the surface 
of last scattering had a perfectly uniform temperature |30| . In particular, the 
fluctuation of the temperature (at the lowest order) in a particular point of the 
space is 

^(6) = i/ d\^-h rr . (4) 

1 * Jnullgeodesic or l 

The integral in Eq. (TJJ is taken over a path of null geodesic which leaves the 
current spacetime point heading off in the direction defined by the unit vector 
Cl and going back to the surface of last scattering of the CMBR. 

Here A is a particular choice of the affine parameter along the null geodesic. 
By using conform coordinates, we have for the metric perturbation 

Sg a b = R 2 (v)h a b (5) 

and r in Eq. ([JJ is a radial spatial coordinate which goes outwards from 
the current spacetime point. The effect of a long wavelength GW is to shift the 
temperature distribution of CMBR from perfect isotropy. As the fluctuations 
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are very small (< AT/T >< 5 x 10" 5 [22K23]), tlle perturbations caused by 
the relic SBGWs cannot be too large. 

The WMAP results give rather tight constraints on the spectrum of the 
SBGWs at very long wavelengths. In [231 Hi] we nn d a constraint on tt gw (f) 
which derives from the COBE observational limits, given by 

*W/)>»ioo < 7 x ^"(y) 2 for H Q < f < 30H Q . (6) 

We included a dimensionless factor h W Q which now is just a convenience (in 
the past it came from an uncertainty in the value of Ho). From the WMAP 
data we have h wo = 0.72 ± 0.05 [2211215]. 

Now, the same constraint will be obtained from the WMAP data [20]. Be- 
cause of its specific polarization properties, relic SBGWs should generate par- 
ticular polarization patterns of CMBR anisotropies, but the detection of CMBR 
polarizations is not fulfilled today [31j . Thus, an indirect method will be used. 
We know that relic GWs have very long wavelengths of Hubble radius size, 
then the CMBR power spectrum from them should manifest fast decrease at 
smaller scales (high multipole moments). But we also know that scalar modes 
produce a rich CMBR power spectrum at large multipole moments (series of 
acoustic peaks, [221 123] )■ Then, the properties of tensor modes in the cosmolog- 
ical perturbations of spacetime metric can be extract from observational data 
using angular CMBR power spectrum combined with large scale structure of 
the Universe. One can see (Fig. 1 ) that in the range 2 < I < 30 (the same used 
in [25], but with the old COBE data [32]) scalar and tensor contributions are 
comparable. 

From [521 [23], the WMAP data give for the tensor/scalar ratio r the con- 
straint r < 0.9. (r < 0.5 in the COBE data [22]); Novosyadly and Apunevych 
obtained Q S caiar(Ho) < 2.7 x 10 -9 [31]. Thus, if one remembers that, at order 
of Hubble radius, the tensor spectral index is —4 < n t < —2 [20], it results 

*W/)>»ioo < L6 x 10" 9 (y) 2 for H <f< 30H , (7) 

which is greater than the COBE data result of eq. ([6]). 

We emphasize that the limit of Eq. ([7]) is not a limit for any GWs, but only 
for relic ones of cosmological origin, which were present at the era of the CMBR 
decoupling. Also, the same limit only applies over very long wavelengths (i.e. 
very low frequencies) and it is far below frequencies of the Virgo - LIGO (Laser 
Intcrfcrometric Gravitational Observatory) band. 

The primordial production of the relic SBGWs has been analysed in [25l 
[26l [27] [33] (note: a generalization for f(R) theories of gravity will be given 
in the next section following |21|). where it has been shown that in the range 
10~ 15 Hz < f < 10 10 Hz the spectrum is flat and proportional to the ratio 

« 10- 12 . (8) 

PPlanck 
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Spherical Harmonics I 

Figure 1: the tensor to scalar ratio, adapted from C. Corda Mod. 
Phys. Lett. A 22, 16, 1167-1173 (2007) 
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Figure 2: the WMAP bound, adapted from C. Corda - Astropart. 
Phys., 30, 209-215 (2008) 



As the spectrum remains flat at high frequency, the relic SBGWs are quite 
important also in the framework of HFGWs, which have been emphasized in 
some works in the recent literature [48] [49J [50] [51] . 

WMAP observations put strongly severe restrictions on the spectrum, as 
we discussed above. In Fig. 2 the spectrum Q gw is mapped: the amplitude 



(determined by the ratio 



PPla 



-) has been chosen to be as large as possible. 



consistent with the WMAP constraint (J7J). 

Nevertheless, because the spectrum falls off oc f~ 2 at low frequencies [25] 
[25] [27] [33J , this means that today, at Virgo and LISA frequencies, indicated in 
Fig. 2, 



n gw (f)hloo <9x 1(T 13 , (9) 
while using the COBE data it was 

ttgw(f)hl 00 < 8 x l(T 14 (refs. [23132]). 
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It is interesting to calculate the correspondent strain at « 100Hz, where 
interferometers like Virgo and LIGO have a maximum in sensitivity. The well 
known equation for the characteristic amplitude [201 [HI [H] can be used: 



h c (f) ~ i.26 x icr 18 (— A)^hj 00 n gw (f), (io) 



obtaining 



h c (W0Hz) < 1.7 x 10" 26 . (11) 

Then, as for ground-based interferometers a sensitivity of the order of 10 -22 
is expected at ss WOHz. four order of magnitude have to be gained in the signal 
to noise ratio [01 EH EQ HI EH Ell E3 HE] ■ Let us analyze smaller frequencies 
too. The sensitivity of the Virgo interferometer is of the order of 10 -21 at 
« 10Hz [HI HHj and in that case it is 

h c (lOHz) < 1.7 x 10~ 25 . (12) 

For a better understanding of the difficulties on the detection of the SBGWs, 
a lower bound for the integration time of a potential detection with advanced 
LIGO is released. For a cross-correlation between two interferometers the signal 
to noise ratio (SNR) increases as 




where Pj(|/|) is the one-sided power spectral density of the i detector [34] 
and 7(/) the well known overlap-reduction function [Ml [35]. Assuming two 
coincident co-aligned detectors (7(/) = 1) with a noise of order 10~ 48 / Hz (i.e. 
a typical value for the advanced LIGO sensitivity [35]), one gets (SNR) ~ 1 for 
~ 3* xlQ 5 years using our result Q gw (f)hl m ~9x 10~ 13 while it is (SNR) ~ 1 
for ~ 3 x 10 7 years using previous COBE result Vt gw (f)h\ m ~ 8x 10~ 14 . As the 
overlap reduction function degrades the SNR, these results can be considered 
a solid upper limit for the advanced LIGO configuration for the two different 
values of the spectrum. 

The sensitivity of the LISA interferometer will be of the order of 10 -22 at 
10~ 3 iJz [37] and in that case it is 

h c (10~ 3 Hz) < 1.7 x 10~ 21 . (14) 

Then, a stochastic background of relic gravitational waves could be in prin- 
ciple detected by the LISA interferometer. We also hope in a further growth in 
the sensitivity of advanced projects. 

We emphasize that the assumption that all the tensor perturbation in the 
Universe are due to a SBGWs is quite strong, but our results ©, (fTTj) , (fT2"]) and 
([T4)) can be considered like upper bounds. 

Resuming, in this section the SBGWs has been analysed with the aid of the 
WMAP data, while previous works in the literature used the old COBE data, 
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seeing that the predicted signal for these relic GWs is very weak. From our 
analysis it resulted that the WMAP bound on the energy spectrum and on the 
characteristic amplitude of the SBGWs are greater than the COBE ones, but 
they are also far below frequencies of the earth-based antennas band. In fact, the 
integration time of a potential detection with advanced interferometers is very 
long, thus, for a possible detection, we have to hope in a further growth in the 
sensitivity of advanced ground based projects and in the LISA interferometer. 



3 Stochastic background and f(R) theories of 
gravity 

GWs are the perturbations h a p of the metric g a p which transform as three- 
tensors. Following [5TJ[35], the GW-equations in the TT gauge are 

Dhj = 0, (15) 

where □ = {—g)~ 1 / 2 d a {—g) 1 l 2 g al3 dp is the usual d'Alembert operator and 
these equations are derived from the Einstein field equations deduced from the 
Einstein-Hilbert Lagrangian density L = R [6l [21] . Clearly, matter perturba- 
tions do not appear in (fT5"|) since scalar and tensor perturbations do not couple 
with tensor perturbations in Einstein equations. The Latin indices run from 1 
to 3, the Greek ones from to 3. Our task is now to derive the analogous of 
Eqs. (|15p assuming a generic theory of gravity given by the action 



1 

21- 



d 4 xV=g~f(R), (16) 



where, for a sake of simplicity, we have discarded matter contributions. A 
conform analysis will help in this goal. In fact, assuming the conform transfor- 
mation 

g a [i = e 2 ®g a /3 (17) 

where the conform rescaling 

e 2 * = \f'(R)\ (18) 

has been chosen being the prime the derivative in respect to the Ricci curva- 
ture scalar and <E> the "conform scalar field" , we obtain the conform equivalent 
Hilbcrt-Einstcin action 



where 3> ;Q ) is the conform scalar field contribution derived from 

Rap = R a p + 2($ ;Q $;/3 - g a ^-s^' s - \g a ^' s -,s) (20) 
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and 

R = e" 2 * + (R - 6D$ - 6$ : i$ ;<5 ). (21) 

In any case, as we will see, the £($, $ ;Q )-term does not affect the GWs-tensor 
equations, thus it will not be considered any longer (note: a scalar component in 
GWs is often considered (TTJ 39, 40, 41j[52], but here we are taking into account 
only the genuine tensor part of stochastic background) . 

Starting from the action (TIT))) and deriving the Einstein-like conform equa- 
tions, the GWs equations are 

Ohl = 0, (22) 

expressed in the conform metric g a p- As no scalar perturbation couples to 
the tensor part of gravitational waves, it is 

h{ = 9 l] Sg u = e- 2 V J e 2 %^ = K, (23) 
which means that h\ is a conform invariant. 

As a consequence, the plane wave amplitude h\ = h(t)el exjp(ikiX J ), where 
e\ is the polarization tensor, are the same in both metrics. In any case the 
d'Alembert operator transforms as 

□ = e- 2 *(D + 2$ ;Q 5 ;Q ) (24) 

and this means that the background is changing while the tensor wave am- 
plitude is fixed. 

In order to study the cosmological stochastic background, the operator (|24[) 
can be specified for a Friedman- Robertson- Walker metric [26j [27] , and Eq. (|22[) 
becomes 

h + (3H + 2$)h + k 2 a- 2 h = 0, (25) 

being □ = + 3if Jj, a(t) the scale factor and k the wave number. It 
has to be emphasized that Eq. (|25|) applies to any f(R) theory whose conform 
transformation can be defined as e 2 * = f'(R). The solution, i.e. the GW 
amplitude, depends on the specific cosmological background (i.e. a(t)) and the 
specific theory of gravity. For example, assuming power law behaviors for a(t) 
and $(t) = |Ln /'(#(*)), that is 

m) = f\R) = f Q {^T\ a(t) = a (h n , (26) 
it is easy to show that general relativity is recovered for m = while 

m 2 + m - 2 

n = — 5— 27 

m + 2 

is the relation between the parameters for a generic f(R) = foR s where 
s = 1 — y with s ^ 1 ,42]. Equation (|25p can be recast in the form 

h+ (3n + m)r 1 /i + fc 2 a ( — ) 2 ™/i = 0, (28) 
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whose general solution is 



h(t) = CjT^CiMx) + C 2 J- a {x)]. (29) 

J a 's are Bessel functions and 

1 — 3n — ui 1 — 3n — m fci 1- " 
a = 77? 7T> P = 7. > x =l > ( 30 ) 

while to, C\, and C2 are constants related to the specific values of n and m. 

The time units are in terms of the Hubble radius H _1 ; n = 1/2 is a radiation- 
like evolution; n = 2/3 is a dust-like evolution, n — 2 labels power-law infla- 
tionary phases and n = —5 is a pole-like inflation. From Eq. (|27[) . a singular 
case is for m = —2 and s = 2. It is clear that the conform invariant plane- wave 
amplitude evolution of the tensor GW strictly depends on the background. 

Now, let us take into account the issue of the production of relic GWs con- 
tributing to the stochastic background. Several mechanisms can be considered 
as cosmological populations of astrophysical sources [43], vacuum fluctuations, 
phase transitions |24| and so on. In principle, we could seek for contributions 
due to every high-energy physical process in the early phases of the Universe 
evolution. 

It is important to distinguish processes coming from transitions like inflation, 
where the Hubble flow emerges in the radiation dominated phase and processes, 
like the early star formation rates, where the production takes place during the 
dust dominated era. In the first case, stochastic GWs background is strictly re- 
lated to the cosmological model. This is the case we are considering here which 
is, furthermore, also connected to the specific theory of gravity. In particular, 
one can assume that the main contribution to the stochastic background comes 
from the amplification of vacuum fluctuations at the transition between an in- 
flationary phase and the radiation-dominated era. However, in any inflationary 
model, we can assume that the relic GWs generated as zero-point fluctuations 
during the inflation undergoes adiabatically damped oscillations (~ 1/a) until 
they reach the Hubble radius H" 1 . This is the particle horizon for the growth of 
perturbations. On the other hand, any other previous fluctuation is smoothed 
away by the inflationary expansion. The GWs freeze out for a/k 3> H^ 1 and re- 
enter the H^ 1 radius after the reheating in the Friedman era (211 [25l [26l [27l [33] . 
Re-entering in the radiation dominated or in the dust-dominated era depends on 
the scale of the GW. After the re-entering, GWs can be detected by their Sachs- 
Wolfe effect on the temperature anisotropy at the decoupling [3D] . When $ 
acts as the inflaton J5TJ 23] we have $ <C H during the inflation. Considering 
the conformal time drj = dt/a, Eq. ([25]) reads 

h" + 2—ti + k' 2 h = 0, (31) 
7 

where 7 = ae* and derivation is with respect to r\. Inflation means that 

a(t) = aoexp(Ht) and then 77 = J dt/a = l/(aH) and ^ = — ~. The exact 
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solution of (|31[) is 



h(r)) = k~ 3/2 y/2/k[Ci(sinkri - coskrj) + C 2 (sinfcr? + cos krj)]. (32) 

Inside the l/H radius it is krj 3> 1. Furthermore, considering the absence of 
gravitons in the initial vacuum state, we have only negative-frequency modes 
and then the adiabatic behavior is 

h = fc 1/2 V^-77 C ' ex PH fc7 7)- (33) 
an 

At the first horizon crossing (aH = k), the averaged amplitude Ah = 
(k/2n) 3 ^ 2 \h\ of the perturbations is 

A h = ±,C. (34) 

When the scale a/k grows larger than the Hubble radius l/H, the growing 
mode of evolution is constant, that it is frozen. This situation corresponds to 
the limit — krj < 1 in Eq. (|32p. Since $ acts as the inflaton field, it is <& ~ 
at re-entering (after the end of inflation). Then, the amplitude Ah of the wave 
is preserved until the second horizon crossing after which it can be observed, 
in principle, as an anisotropy perturbation of the CBR. It can be shown that 
~ A h is an upper limit to Ah since other effects can contribute to the 
background anisotropy [45] . From this consideration, it is clear that the only 
relevant quantity is the initial amplitude C in equation (|33[) which is conserved 
until the re-entering. Such an amplitude directly depends on the fundamental 
mechanism generating perturbations. Inflation gives rise to processes capable of 
producing perturbations as zero-point energy fluctuations. Such a mechanism 

/XT 

depends on the adopted theory of gravitation and then -=r- could constitute a 
further constraint to select a suitable f(R)-thcory. Considering a single graviton 
in the form of a monochromatic wave, its zero-point amplitude is derived through 
the commutation relations 

[h(t,x),n h (t,y)]=iS 3 (x-y), (35) 

calculated at a fixed time t, where the amplitude h is the field and 7T/j is the 
conjugate momentum operator. Writing the Lagrangian for h 

L = \^g aP h- a h,p (36) 
in the conformal FRW metric g a p (h is a conform invariant), we obtain 

n h = g = e 2 *a 3 h. (37) 
oh 



Equation ([55]) becomes 



[h(t,x)My,y)] = i S3 %J ) (38) 
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and the fields h and h can be expanded in terms of creation and annihilation 
operators 

h(t,x) = J d 3 k[h(t)e- lkx + h*(t)e ik % (39) 

h(t, x) = 1 / d 3 k[h(t)e- lkx + h*{t)e lkx ]. (40) 



(2tt) 3 / 2 

The commutation relations in conform time arc then 

[ h { h '*-h*h']=i^. (41) 

Putting ((33]) and ([33]). we obtain C = \/27r 2 J ffe _ * where iJ and $ are 
calculated at the first horizon crossing and then 

A h = ^He~*, (42) 

which means that the amplitude of GWs produced during inflation directly 
depends on the given f(R) theory being $ = | In \f'(R)\. Explicitly, it is 

A h = ; H (43) 
VWW\ 

This result deserves some discussion and can be read in two ways. From one 
side the amplitude of relic GWs produced during inflation depends on the given 
theory of gravity that, if different from general relativity, gives extra degrees of 
freedom which assume the role of inflaton field in the cosmological dynamics [44] . 
On the other hand, the Sachs- Wolfe effect related to the CMBR temperature 
anisotropy could constitute a powerful tool to test the true theory of gravity at 
early epochs, i.e. at very high redshift. This probe, related with data a medium 
[46] and low redshift [47] , could strongly contribute to reconstruct cosmological 
dynamics at every scale, to further test general relativity or to rule out it against 
alternative theories and to give constrains on the SBGWs, if f(R) theories ares 
independently probed at other scales. 

For a better understanding of the issue, it is important to analyse in spirit 
of our discussion some particular viable model ,53], for example the interesting 
viable models in [5j I55| 156] . In this way we will also integrate an analysis which 
started in [57] . 

Let us start by the model proposed in [56] [57] . In this case it is 

Jy ' l + 1 R n ' v ' 

which, in our case, gives the conform scalar field 



2 1 (l + 7 i?"^ 2 
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Here a, (3, 7 are positive constants and n is a positive integer [5BJ [57] . These 
equations cannot satisfy the curvature Rq in the present Universe for which it is 
/(i?o) = because this implies a singularity in the amplitude of the relic GWs 

dlSD. 

In the model proposed in [5] it is 

f(R) = -a[tanh b(R - Ra) + tanh (46) 

and 

$ = -ln| — 1 (47) 

2 ln| 2cosh 2 ^il l [ ' 

Here a, and b are positive constants [51 [57]. Thus, in the case of flat Universe 
R = 0, it is /(0) = and 



cosh 2 M*_iM 
A h = H] j\ _^_|. (48) 

On the other hand, in the case in which it is R Ss> Rq in the present Universe 
it is [5J [57] 

/(#)-► 2A e// = -a[H-tanh^], (49) 

i.e. A e // plays the role of an effective cosmological constant. In this case 
Ah increases with increasing curvature. 
Ah has a maximum which is 



A h =H^\ (50) 
in the case in which R = Rq. 

Two pretty realistic f (R) theories have been proposed in [55J . The first is 



(R-R ) 2n+1 +R 2 n+1 
fo + fi{(R-Ro) 2n+1 +R 2 n+i y 



which gives the correspondent conform scalar field 

2 l (/o + /i{( J R-i?o) 2 " +1 +^" +1 }) 2 ' V ' 

Here n is a positive integer fi = 1/Aj, where A, represents the effective 
cosmological constant and /o ~ 5-R 2 ™ [55J. In this case the amplitude of the 
stochastic background admits a minimum at the point R = R, where R satisfies 

£ I £ p 2n+l 
(5 n \2n+l JO T jynp , 

{R - Ro) - {n+l)h (53) 
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and the minimum is given by 



2 

The second proposal in [55J is the theory 

R 



Ab , H ,hJl^±. ( 64 ) 



//? ^ n n 

dR ^~ a w^r^w^ (55) 



which gives the correspondent conform scalar field 

1 D 2n n -i 

^2^(iTlkF"/^ )+ 2 ln/0 - (56) 

Here a, (3, fo and i?iare constants, while /(+oo) ~ fofiAi is identified with 
the effective cosmological constant at the inflation era if /o/3 = 1. The case of 
flat Universe is interesting. We get 



*~^|=^p|. (57) 

Resuming, in this section it has been shown that amplitudes of tensor GWs 
are conform invariant and their evolution depends ond the cosmological SBGWs. 
Such a background is tuned by a conform scalar field which is not present in 
the standard general relativity. Assuming that primordial vacuum fluctuations 
produce a SBGWs, beside scalar perturbations, kinetic distortions and so on, 
the initial amplitude of these ones is function of the f(R)-theory of gravity and 
then the SBGWs can be, in a certain sense, "tuned" by the theory. Vice versa, 
data coming from the Sachs- Wolfe effect could contribute to select a suitable 
f(R)-theory which can be consistently matched with other observations. 

We have also considered some interesting examples which have been recently 
discussed in the literature [5J [55J 15S1 HI] . 

However, further and accurate studies are needed in order to test the relation 
between Sachs- Wolfe effect and f(R) gravity. This goal could be achieved in the 
next future through the forthcoming space (LISA) and ground based (Virgo, 
Advanced L1GO) interferometers and Li-Baker HFGW detectors. 



4 Conclusions 

This paper has been a review of previous work on the SBGWs which has been 
discussed in various peer-reviewed journals and international conferences. The 
SBGWs has been analysed with the aid of the WMAP data while, in general, in 
previous works in the literature about the SBGWs, old COBE data were used. 
After this, it has been shown how the SBGWs and f(R) gravity can be related, 
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showing, vice versa, that a revealed SBGWs could be a powerful probe for a 
given theory of gravity. In this way, it has also been shown that the conform 
treatment of SBGWs can be used to parametrize in a natural way f(R) theories. 
Some interesting examples which have been recently discussed in the literature, 
see [SJ [SH [SU [57] , have been analysed. 

The presence and the potential detection of the SBGWs is also quite im- 
portant in the framework of the debate on high-frequency gravitational waves 
(HFGWs). The importance of HFGWs has been recently emphasized in some 
papers in the literature and the ground based Li-Baker HFGW has been pro- 
posed to be sensitive at 10 10 GHz to 1CT 32 amplitude HFGWs gS]. 
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